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The trigonal zigzag nanodisk with size N has TV localized spins. We investigate its thermodynam- 
ical properties with and without external leads. Leads are made of zigzag graphene nanoribbons or 
ordinary metallic wires. There exists a quasi-phase transition between the quasi-ferromagnet and 
quasi-paramagnet states, as signaled by a sharp peak in the specific heat and in the susceptability. 
Lead effects are described by the many-spin Kondo Hamiltonian. A new peak emerges in the spe- 
cific heat. Furthermore, the band width of free electrons in metallic leads becomes narrower. By 
investigating the spin-spin correlation it is argued that free electrons in the lead form spin-singlets 
with electrons in the nanodisk. They are indications of many-spin Kondo effects. 
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Introduction: Graphene nanostructurepj has the po- 
tential for future application in nanoelectronics and spin- 
tronics. In particular, much attention is now focused 
on graphene nanoribbons 0] due to almost flat low- 
energy band at the Fermi level depending on the edge 
states. Another basic element of graphene derivatives is 
a graphene nanodiskj^, H, H, 0,1^ ■ It is a nanometer- 
scale disk-like material which has a closed edge. There 
are many type of nanodisks, among which the trigonal 
zigzag nanodisk is prominent in its electronic and mag- 
netic properties because there exist iV-fold degenerated 
half-filled zero-energy states when its size is N (Fig[T]). 

In this paper we explore thermodynamical properties 
of the trigonal zigzag nanodisk. The system is well ap- 
proximated by the infinite-range Heisenberg model. It 
is exactly solvable. A sharp peak emerges in the spe- 
cific heat and in the susceptibility, which we interpret as 
a quasi-phase transition between the quasi-ferromagnet 
and quasi-paramagnet states. We then investigate a 
nanodisk-lead system, where the lead is made of a zigzag 
graphene nanoribbon or an ordinary metallic wire. We 
refer to it as a graphene lead or a metallic lead. 

It is shown that lead effects are described by the many- 
spin Kondo Hamiltonian. Electron spins in the nanodisk 
and the lead orient into the opposite directions to lower 
the coupling energy. A new peak appears around a cer- 
tain temperature (T — Tk) in the specific heat but not 
in the susceptibility for small size nanodisks. The inter- 
nal energy decreases near the peak, and the band width 
of free electrons in the lead becomes narrower in the in- 
stance of the metallic lead. Furthermore, the spin-spin 
correlation takes the maximum value at T = 0, remains 
almost constant for T < Tk, and then decreases mono- 
tonically as T increases. We interpret these phenomena 
to mean that free electrons in the lead are consumed to 
make spin-singlets with electrons in the nanodisk. It is 
intriguing that all electrons (only a few portion of elec- 
trons) in the nanodisk are engaged in the case of the 
graphene (metallic) lead. They are indications of Kondo 
effects due to the Kondo coupling between electrons in 
the lead and the nanodisk. 
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FIG. 1: (a) Trigonal zigzag nanodisks. The size parameter A'^ 
is defined in this way. The number of carbon atoms is given 
by Nc = + &N + 6. (b) The nanodisk-lead system. The 
nanodisk with A = 7 is connected to the right and left leads. 



Quasi- Ferromagnet: The size-A^ zigzag trigonal nan- 
odisk has A^-fold degenerated zero-energy states [4], where 
the gap energy is as large as a few eV. Hence it is a good 
approximation to investigate the electron-electron inter- 
action physics only in the zero-energy sector, by project- 
ing the system to the subspace made of those zero-energy 
states. The zero-energy sector consists of N orthonormal 
states \fa), a = 1,2,--- ,N, together with the SU(N) 
symmetry. Let Uaf3 and Jap be the Coulomb energy and 
the exchange energy between electrons in the states \fa) 
and If,}). It follows (yj that Jaf3 ~ Uafi and that all Jaf3 
are of the same order of magnitude for any pair of a and 
/3, implying that the SU(N) symmetry is broken but not 
so strongly. It is a good approximation to start with the 
exact SU(A^) symmetry, by replacing Uafi and Jap with 
their averages U and J, respectively. Then, the zero- 
energy sector is described by the Hamiltonian [9], 



Hy 



-jsi 



+ J ntot, 



(1) 



with [/' = [/ ~ i J and J « U, where Stot = Ea ^ («) 
is the total spin, and n^ot — i'^) ^he total elec- 

tron number. Here, 71(0;) = dl.{a)d„{a), and S{a) — 
^dl.{a)Ta-a'da'{oi) with da{a) the annihilation operator 
of electron with spin a =t, J. in the state \fa)- t is the 
Pauli matrix. We call Stot the nanodisk spin. 

Apart from an irrelevant constant, the Hamiltonian ([1]) 
is the infinite-range Heisenberg model, Hs = — JSjot: in 
the half-filled sector with ntot — N. The nanodisk spin 
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FIG. 2: Thermodynamical properties of the nanodisk without 
leads, (a) The specific heat C in unit of k^N . (b) The entropy 
S in unit of fcBA''log2. (c) The magnetization (Sjot) in unit 
of Sg. (d) The susceptibihty x in unit of Sg. The size is A'^ = 
1, 2, 2^, • • ■ 2^*^. The horizontal axis stands for the temperature 
T in unit of JN/ks ■ The arrow represents the phase transition 
point Tc in the limit TV — > oo. 



takes the maximum value Sg — \J N/2 {N/2 + 1) in the 
ground state, where all spins are spontaneously polarized. 
The nanodisk spin system exhibits a strong ferromagnetic 
order due to a large exchange interaction. The relaxation 
time is finite but quite large even if the size N is small. 
We have called such a system quasi-ferromagnet[l|. 

Thermodynamical Properties: The infinite-range 
Heisenberg model is exactly diagonahzable, iJsl'^'s) — 
Es\'i!s), with Es = -JS{S + 1), where S takes half- 
integer or integer values from N/2 down to 1/2 or 0, 
depending on whether N is odd or even. The total de- 
generacy of the energy level Es is (25*+ l)gfq{S) with 

gN {N/2 - q) =NCq -NCq-l. 

We have a complete set of the eigenenergies together 
with their degeneracies. The partition function of the 
nanodisk with size N is exactly calculable, 



Z^=Y,{2S+l)g^{S)e-^'^^^+'\ 



(2) 



According to the standard procedure we can evaluate the 
specific heat Cn(T'), the entropy S^{T), the magneti- 
zation (S?ot) and the susceptibility x = fe^((S?ot) " 
(Stot)'^) from this partition function. The entropy is 
given by 5n (0) = fee log(7V -|- 1) at zero temperature, 
corresponding to the ground state multiplicity -I- 1. 
We display them in FiglUfor size iV = 1, 2, 2^, • • • 2^°. 

There appear singularities in thermodynamical quan- 
tities as ^ oo, which represent a phase transition at 
Tc = JN/2kB between the ferromagnet and paramag- 
net states (Figl2|). For finite N, there are steep changes 
around Tc, though they are not singularities. It is not a 
phase transition. However, it would be reasonable to call 
it a quasi-phase transition between the quasi-ferromagnet 



and quasi-paramagnet states. Such a quasi-phase transi- 
tion is manifest even in finite systems with N ~ 100 ~ 
1000. 

The specific heat and the magnetization take nonzero- 
values for T > Tc [Fig[2{a),(c)], which is zero in the limit 
N ^ oo. The entropy for T > Tc is lower than that of the 
paramagnet [Figllt^b)] . These results mean the existence 
of some correlations in the quasi-paramagnet state. The 
maximum value of the susceptibility increases linearly as 
TV becomes large. It is an indicator of the quasi-phase 
transition. 

Many-Spin Kondo Hamiltonian: We proceed to inves- 
tigate how thermodynamical properties of the nanodisk is 
affected by the attachment of the leads. Though there are 
two leads attached to a nanodisk, the lead Hamiltonian 
i?L and the transfer Hamiltonian Ht are expressed as if 
there were a single lead after a certain transformation^], 

Hi^ = J2^('^)4aCka, (3a) 

= iJ2 E {4.dAc^) + 4(«)cfc.) , (3b) 

where c^a- is the annihilation operator of electron in the 
lead with the wave number k and the dispersion relation 
e{k). 

When charges transfer between the nanodisk and the 
leads, the total electron number ntot is no longer fixed 
in the nanodisk. However, the nanodisk remains to be 
half-filled, when a charge transfers from the lead to the 
nanodisk and then transfers back from the nanodisk to 
the lead. The process is the second order effect in the 
tunneling coupling constant t. We derive the effective 
Hamiltonian for such a process. 

The total Hamiltonian is H = Hd + + Ht. 
We take Hq = Hu + Hi, as the unperturbed term 
and Ht as the perturbation term. Note that U » 
t. We make a canonical transformation known as 
the Schrieffer- Wolff transformationjl^ to eliminate H^, 
H ^ H = e*'^iJe^''^, with G the generator satisfying 
Ht + I [G, Hq] = 0. The dominant contribution comes 
from the Fermi surface, e (fc) = ep- We assume the sym- 
metric condition ep = U + ■^U' with respect to the Fermi 
energy. Then, after a straightforward calculation, we ob- 
tain H = H-o + H\, + Hk + O (P ) , where Hk is the second 
order term in <. It is the many- spin Kondo Hamiltonian, 



EIk = Jk 



E 

kk' (T(t' 



„t 



'Cfc', 



•St 



(4) 



with the Kondo coupling constant Jk = 8P/U' . The 
difference between the above many-spin Kondo Hamil- 
tonian and the ordinary Kondo Hamiltonian is whether 
the local spin is given by the summation over many spins 
Stot or a single spin S. Note that S^^^ is a dynamical 
variable but is not, = 3/4. 
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FIG. 3: Thermodynamical properties of the nanodisk (a,b) 
with graphene leads and (a',b') with metallic leads. See the 
caption of Figlg] We have set Jk/J = 0.2 and D = 2kBTc. A 
new peak appears around Tk in the specific heat. 



Functional Integration: The total Hamiltonian is now 
given by iJoff — Hs + ^^l + Hk at half filling. We de- 
fine the spinor = (c|,C|)*. The partition function in 
the Matsubara form is given in terms of the Hamiltonian 
density TCcS as 



Z = Trs J ViJjVip^ exp 
-Trs [exp(-/3i/s)^K 



J^^ J dx {i^'^dril^ + n^s) 



(5) 



with ^'k = / VipV^^ exp[— S'k], where S'k is the action 
5k = dr J dx {i;''drip + + Hk)- We first perform 
a functional integral over the lead electron's degree of 
freedom in Zk, and then summed up the nanodisk spin 
in (O. 

Because an electron in the lead is constrained within a 
very narrow region, it is a good approximation to neglect 
momentum scatterings, 



Hk — 2 Jks • St 



(6) 



where s — ^ J2kaa' cj.o.''o-<T'Cfc<j' is the electron spin in the 
lead. The action S'k is summarized as 



duj 
2^ 



(fc) M (k) ^ (k) 



(7) 



with M (k) — —[iuj ~ £ (k)] + Jkt ■ ^tot- Perform- 
ing the functional integration we find Zk =Det[M] = 
exp [—/3Fk], where Fk is the Helmholtz free energy, 

FK^~^J2^n [cosh(/3JK IStotI) + cosh(/3e (fc))] . (8) 



This is reduced to the well-known formula for free elec- 
trons with the dispersion relation e(fc) for Jk = 0. 



Zigzag Graphene Nanoribbon Leads: We first consider 
the system where the leads are made of zigzag graphene 
nanoribbons. Owing to the flat band at the zero energy, 
e (fc) = 0, the result of the functional integration ([8]) is 
quite simple. 



^^K = In cosh[/3JK IStotI /2]. 



(9) 



The effective Hamiltonian for the nanodisk spin is Hs + 
Fk- The lead effect is to make the effective spin stiffness 
larger and the ferromagnet more rigid. 

The trace over the nanodisk spin is carried out in ([5|). 

Z = ^ (25 + 1) gjv(5)e"^^^(^+i) ^^^^ (^^^/s{S + T 
s 

(10) 

We compare thermodynamical properties of the nanodisk 
with leads (Fig[3]) and without leads (Fig[2]). The mag- 
netization (S^Qt) and the susceptibility x are found to 
be indistinguishable from those of the nanodisk without 
leads (Figl21). In FigO we show the specific heat Cg{T) 
and the entropy Sg{T) for various size N. The significant 
feature is the appearance of a new peak in the specific 
heat at T = Tk « (Jk/2J)Tc, though it disappears for 
large N. We examine the internal energy Eg{T), which 
is found to decrease around Tk (FigH]). Near zero tem- 
perature it reads 



EciT) -JSl - ^5g 



JkSqC 



(11) 



The first term (oc J) represents the energy stabilization 
due to the ferromagnetic order present in the nanodisk 
without leads, while the second term (oc Jk) represents 
the one due to the Kondo coupling. Furthermore, it fol- 
lows that the entropy is reduced at zero temperature as 
Sg (O) — Sn (0) = —fee log 2, as imphes that the ground 
state multiplicity at the zero temperature is just one half 
of that of the system without leads. 

The spin-spin correlation |(s • Stot)| is calculated based 
on the partition function (O and shown in FiglSja). Near 
zero temperature we find 



1 



|(s-Stot)| ^ -'5'gtanh [/?JKSg/2] 



(12) 



It takes the maximum value S'g/2 at T = 0, and remains 
almost constant for T < Tk, and then monotonically 
decreases as T increases. Finally, it almost vanishes in 
the quasi-paramagnet phase for large N since (Stot) ~ 0. 
We may interpret these phenomena as follows. Electrons 
in the lead and the nanodisk form spin-singlet states to 
lower the coupling energy ^ . The singlet state is rather 
tight for T <Tk, but thermally broken as T increases. 

All these features indicate the occurrence of the Kondo 
effect due to the coupling between the spins in the nan- 
odisk and the leads. 
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(a) without leads 

0.2 0.4 +"^0.6 



(b) with graphene leads 




FIG. 4: The internal energy E in unit of Nku for the nan- 
odisk (a) without leads and (b) with graphene leads. See the 
caption of Fig[2] (a) The energy descreases except for A'^ = 1 
as the temperature decreases, which represents the ferromag- 
netic order, (b) There exists an additional energy decrease 
around Tk, which is prominent for A'' = 1, attributed to the 
Kondo effect. 



(a) with graphene leads 



(b) with metallic leads 




FIG. 5: The spin-spin correlation in unit of Sg/2. See the 
caption of Fig|3] The correlation occurs due to the Kondo 
coupling. It is almost constant for T < Tk, and decreases as 
T increases. 



Metallic Leads: Next we consider the system com- 
prised of metallic leads with a constant energy density, 
p(e) = p for |e| < D and p{e) — Q for \e\ > D. We 
change the momentum integration into the energy inte- 
gration in ([8]), 



"2/3 



I deln [cosh/JJx |Stot| + cosh/3e] . (13) 

J-D 



The free energy is given by Fk = -p/3~^ |Stot| - 

2/3Dlog2-hLi2( - e''(^"'^Kl^'°'l))-Li2( 
g-/3(D-i-JK|Stot|)^j ^ where Li2 [x] is the dilogarithm 
functionfm. It is reduced to the free energy of the 
nanodisk with graphene leads in the limit £) — > with 
p= 1/2D. 

The magnetization (StQ^.) and the susceptibility x ^-re 
found to be indistinguishable from those in the case of no 
leads (FigUl). The relation Sm (0) - S'n (0) = -fee log 2 
holds for the entropy precisely as in the case of graphene 
leads. In FiglH we show the specific heat Cm(T) and the 
entropy Sm{T) for various size A^. A broad peak appears 
at T = Tk in the specific heat. Using the asymptotic 
behaviors [Tlj of Li2 [a;], we obtain the free energy in low 



temperature regime as 

/3i?JK|Stot|- 



7r2 p , 



Jk I St 



(14) 

which is valid up to the terms in the order of e . The 
entropy Sm (T) , the specific heat Cm {T) and the internal 
energy Em (T) read as follows. 



(T) : 

and Em (T) 
AE{T) 



kB log 



iV+ 1 TT^ 



pklT, 



Cm{T) 



-pk 



Eg (T) + AE{T) with 



-^{D-JKS,f + ^p{kBTf. 



2 J, 



(15) 



(16) 



All terms proportional to p have arisen from free elec- 
trons in the metallic lead. The internal energy Em (T) 
consists of two terms: Eq (T) is identical to the energy 
(fTTj) for the nanodisk with graphene leads, and AE{T) is 
the energy of the metallic lead. The first term of AE{T) 
shows that the band width of free electrons in the lead 
becomes narrower due to the Kondo coupling. We may 
interpret that n free electrons with 



pJkSq 



(17) 



are consumed to make spin-coupling with electrons in the 
nanodisk. The second term is the thermal energy of free 
electrons in the metallic lead. 

We show the spin-spin correlation | (s- Stot) | in FiglSl^b) . 
The overall features are the same as for the nanodisk with 
graphene leads. However, there are new features. First 
of all, the value of correlation is quite small. This is 
because spin-singlets are formed only by a small portion 
of electrons in the metallic lead which are near the Fermi 
level. We expect that this number density is given by 
PT)) . Indeed, it is observed that |(s • St,ot)| ^ ^pJkSq at 
T = as N ^ oo [FigEtb)]. 

In this paper we have investigated thermodynamical 
properties of a zigzag graphene nanodisk without and 
with leads. The lead effects are summarized by the 
many-spin Kondo Hamiltonian. One effect is to enhance 
the ferromagnetic order. This result is important to 
manufacture spintronic circuits by connecting leads in 
nanodevicesjil. We have shown various thermodynami- 
cal results indicating many-spin Kondo effects. 
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